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We present a scheme to cool the motional state of neutral atoms confined in sites of an optical 
lattice by immersing the system in a superfluid. The motion of the atoms is damped by the gen- 
eration of excitations in the superfluid, and under appropriate conditions the internal state of the 
atom remains unchanged. This scheme can thus be used to cool atoms used to encode a series of 
entangled qubits non-destructively. Within realisable parameter ranges, the rate of cooling to the 
ground state is found to be sufficiently large to be useful in experiments. 
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I. INTRODUCTION 

Neutral atoms are one of the most promising candi- 
dates as carriers for the storage and manipulation of 
quantum information Q]. Qubits may be stored in long- 
lived internal atomic states with very low levels of de- 
coherence, and may be manipulated using interactions 
between the atoms and external devices (such as lasers) 
or interactions amongst the atoms themselves. 

An experimental prerequisite for this is the develop- 
ment of techniques to trap single atoms, and there has 
been much progress over the last five years both in optical 
traps 0, and in magnetic microtraps . In addition, 
specific implementation of quantum computing usually 
requires cooling of atoms to the vibrational ground state 
of the trap, or at least to the Lamb-Dicke limit. Many 
techniques have been developed including the widespread 
use of laser cooling [j| . 

One of the most promising routes to quantum compu- 
tation with neutral atoms is the use Bose Einstein Con- 
densates (BECs) @ loaded in optical lattices 0, @> a 
system which has been realised in part in a number of 
recent experiments 0, 0, 0, 0] . There are several theo- 
retical proposals for the implementation of quantum logic 
gates in such systems [Tl 111 Hi liH . and the first 
steps towards the fundamental experimental techniques 
required for some of these have been recently realised. 
For example, the recent demonstration of spin-dependent 
transport in an optical lattice [t| makes possible the im- 
plementation of a fundamental quantum phase gate by 
cold controlled collisions 12] in which qubits are encoded 
using two different internal states of the atoms in the op- 
tical lattice. 

However, most of these proposals require the trans- 
port of qubits, which is usually associated with heating 
of the atomic motion 0. A question then arises as to 
how that motion may be cooled back to the ground state 
without changing the internal state of the atoms, and 
thus destroying the qubits or their entanglement. Laser 
cooling, for example, is clearly not applicable here as the 
process of light scattering causes decoherence. The same 
problem arises in scalable ion trap quantum computing, 
and there it has been overcome using sympathetic cooling 
schemes, in which ions used to encode qubits are cooled 



via a coulomb interaction with either a single ion which is 
directly laser-cooled ^7$ or another species of ions which 
are directly laser-cooled ^8|. I n a different context, sym- 
pathetic cooling schemes are also used widely in the field 
of cold quantum gases, where they have been used to 
cool different spin states of the same atomic species , 
to cool different Bosonic species [2j}, and to cool Fermi 
gases brought into contact with a BEC [2l| . 

In this article we consider the sympathetic cooling of 
a single atom in a harmonic trap in contact with a su- 
perfluid. This is readily expanded to the case of many 
harmonic traps, which is a good approximation for an op- 
tical lattice without tunneling. The motion of the atom 
is damped by the generation of excitations in the super- 
fluid, and the resulting cooling rates are sufficiently large 
to be useful experimentally. In addition, decoherence of 
a qubit encoded on the atoms can be eliminated in this 
scheme provided that the internal atomic states used to 
encode the qubit are chosen carefully in order to satisfy 
particular collisional requirements. 



II. OVERVIEW 

In this section we give a short summary of the most 
important results contained in this article. Derivations 
and further discussions of these results follow in the re- 
maining sections. 

Our goal is to cool a single trapped atom represent- 
ing a qubit |0), |1) without destroying the superposition 
state of the qubit (or the entangled state in case of many 
atoms). Cooling of the atom is achieved by sympathetic 
cooling, immersing the atom in a superfluid, which plays 
the role of a very cold reservoir. By a proper encoding 
of the qubit in internal atomic states, and choice of the 
atomic level for the superfluid reservoir (see section llll A|) 
we can ensure that (i) the qubit is not destroyed by open- 
ing collisional channels to unwanted final states, and (ii) 
the |0) and |1) states have identical collisional proper- 
ties with respect to the collisional interactions with the 
superfluid, and thus the collisions do not randomise the 
relative phases of the qubit. 

Cooling is considered within a model in which the 
atoms are treated as being trapped in independent 
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ID Harmonic oscillator potentials with trapping fre- 
quency uj, and interact with the superfiuid via a density- 
density interaction, generating excitations in the super- 
fluid, which are modeled as Bogoliubov excitations in a 
weakly interacting Bose gas (section and have mo- 
mentum ?iq and energy e q . In discussing this cooling 
process we can restrict ourselves to a single component 
of the qubit |0) (or |1)). 

A master equation is derived for the density operator of 
this system C section 1111 CI and appendix [XJ , from which 
the time evolution of the probability p n that the atom 
is in the nth motional state of the Harmonic oscillator 
potential is shown to be 



Pm, — ^ ^ F n — >mPn ^ ^ -^m — >n'Pm~ i r '} ^ Hn.mjPn Pn 



n' <m 



(1) 



The terms with coefficient F n ^ m , 

= ^-J2 l^n,m(q)| 2 <5(Mn - m) - e q ), (2) 



Fn 



where Z„. m are the matrix elements of the interaction 
Hamiltonian in the basis of Harmonic oscillator energy 
eigenstates (Fock states), describe the transitions from 
state n to state m due to generation of excitations in the 
superfiuid, and the terms with coefficient if n>m , 

H m ,n = -r ^iV(q)|Z rMn (q)| 2 (5(Mn-m| - £,), (3) 



describe the transitions between state n and state m due 
to interactions with thermal excitations at finite temper- 
atures. This is illustrated in Fig. 
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FIG. 1: The motion of an atom in a harmonic trap immersed 
in a superfiuid is cooled by the generation of excitations in the 
superfiuid. Terms in the equations of motion with coefficients 
F n ^ m describe transitions from oscillator state \n) to state 
\m) by creation of excitations, whilst finite temperature con- 
tributions with coefficients H n<m account for the interaction 
of the system with thermal excitations. 

If the speed of sound in the superfiuid is u, and the 
mass of superfiuid atoms is nib, then the behaviour of 
the cooling process can be separated into two regimes 
- where the motion of the oscillating atom is subsonic 
(hui <C rribu 2 /2) or supersonic (hcu ^> m^u 2 /2). In the su- 
personic regime, cooling from any excited oscillator state 
occurs directly to all lower energy states, including a sig- 
nificant transition rate directly to the ground state (sec- 
tion IIV A"|) . The resulting rate of energy loss e(n) for a 



particle in the nth oscillator state is not linear in n, but 
instead (for lattice and superfiuid atoms of equal mass 
m) is found to be 



e(n) 



7 3 / 2 



-a[e(n)} 



3/2 



(4) 



where g a b is the coupling constant for interactions be- 
tween the atoms in the lattice and the superfiuid, pa is 
the condensate density, e(n) = hum, and a ~ 0.3 is a 
constant. 

If we consider the slowest transition rate, that from the 
first excited state to the ground state, we find that the 
characteristic transition time, r, is given by 

1 Clnh 

(5) 



2tt 



1.2 x 10" 



Poa 3 ab l ' 



where a a b is the scattering length for the interaction 
between atoms in the lattice and the superfiuid, and 
Iq = yjh/{mus) is the size of the harmonic oscillator 
ground state. Thus, as Iq is typically an order of magni- 
tude larger than a a b, and poa^ ~ 10~ 4 , r is of the order 
of 10 oscillator cycles. This is a sufficiently rapid cooling 
rate to be useful experimentally. 

In the subsonic regime, damping still occurs because 
the oscillatory motion of the atom is accelerated (section 



HVBIl . Significant rates are found only for transitions 
between neighbouring oscillator levels, and e(n) is found 
to be linear in n, 



s(n) 



je(n). 



(6) 



12irm a mbU 

When the superfiuid is at a finite superfiuid tempera- 
ture T, the system is cooled to the temperature T of the 
superfiuid. The final distribution of occupation probabil- 
ities is shown to be a Boltzmann distribution, 

p n = p e~ nhuj/{kBT} = (l - e"^/( fesT ^ e -nHu>/(k B T)^ 

(7) 

where ks is the Boltzmann constant (section llV C|l . If 
the temperature corresponds to an energy much smaller 
than the Harmonic oscillator spacing, ksT fko, then 
the population in excited motional states is negligible. 
For example, if uj ~ 2?t x 10 5 s~ 1 , TiLo/hs ~ 5/iK, so that 
for T = 500nK, we then obtain 1 — po w 5 x 10~ 5 . 

A semiclassical treatment of this system in the WKB 
approximation (section gives a result for the super- 
sonic case which is different from the full quantum result 
by only 12%. A similar treatment in the strongly sub- 
sonic regime gives exact agreement with the earlier result. 

In section I VII we investigate a somewhat different 
model for the excitations, in the context of a quasi-lD 
superfiuid. The resulting damping rates are found to be 
small except in the regime where the superfiuid is very 
strongly interacting, which is a difficult regime to obtain 
experimentally. Finally, appendix [B] contains semiclassi- 
cal estimates for small additional damping terms which 
arise at finite temperatures and have been neglected in 
earlier calculations. 
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III. THE MODEL 
A. Avoiding Decoherence 

The total Hamiltonian of the cooling process for a sin- 
gle atom can be written as 

Hot = ffqubit + -^motion + -^superfluid + Hi nt , (8) 

where -ff qu bit is the Hamiltonian for the internal states 
of the atom, denoted |0) and |1), on which the qubit 
is encoded, -Hmotion is the Hamiltonian for the atomic 
motion of the atom which is to be cooled, -ff SU perfluid is 
the Hamiltonian for the superfluid, and H lnt describes 
the interaction between the atom and the superfluid. In 
order to cool a qubit without decoherence, the internal 
state of the atom being cooled should remain unchanged 
during the cooling process. If we write the initial inter- 
nal state of an atom in a particular lattice site as \ip), 
and the combined density operator for the initial mixed 
motional state of the atom and the state of the super- 
fluid as W-^O), so that the total initial density operator 
is IV'X^/'I <8> W(0), then the overall Hamiltonian for the 
cooling process, //tot, must satisfy 

Thus Htot must be of the form H tot = (|0)(0| + |1)(1|) ® 
H. This requirement is satished provided that the inter- 
action Hamiltonian, H^t is independent of the internal 
state of the atom in the lattice. Thus, the trap poten- 
tial must be the same for the two internal states |0) and 
|1), and the scattering length a a b between atoms in the 
superfluid and atoms in the lattice must also be the 
same for the two internal states. The identical scatter- 
ing lengths can be arranged by choosing symmetric spin 
configurations, for example, by choosing |0) and |1) to be 
internal states with angular momentum quantum num- 
ber F = 1 and magnetic quantum numbers mp = ±1, 
and the superfluid atoms to be in an internal state with 
F = 1 and mp = 0. In order to make such a configura- 
tion stable against spin-exchanging collisions |23J, these 
states should all be in the ground state of the manifold, 
and to prevent the creation of pairs of lattice atoms from 
superfluid atoms j^, the energy of the mp — level 
should be lowered with respect to the mp = ±1 states 
(for example by using a laser j25|). We must also ensure 
that when we have N qubits (N > 1), the entanglement 
between them is not destroyed when the motion of one 
or more of them is cooled. The condition in JjJJl is once 
again sufficient for the suppression of decoherence, but 
now is the total internal state of the iV-qubit sys- 
tem, and W is the total combined density operator for 
the motional state of each qubit and the state of the en- 
vironment. Physically, the condition is now modified so 
that the interaction between any atom and the super- 
fluid must be both independent of the internal state of 
that atom, and independent of the internal state of all 



other atoms. Because the interaction is a density-density 
interaction, this second requirement is always fulfilled. 
Note that when the correlation length of the superfluid 
is shorter than the separation between atoms, it is pos- 
sible for the motional state of different atoms to become 
entangled. However, this will not affect the state of the 
./V-qubit system, as the qubits are encoded solely on the 
internal states of the atoms, which remain at all times 
separable from the motional states. 

B. Hamiltonian for the Oscillator-Superfiuid 
Interaction 

After imposing the requirement from the previous sec- 
tion, we consider only the motional degrees of freedom 
of the atoms in the optical lattice, which are assumed to 
be confined in particular lattice sites where the motional 
states can be approximated as those of an harmonic oscil- 
lator. Coupling to the superfluid occurs in the form of a 
density-density interaction that generates excitations in 
the superfluid, which we model as Bog oliubov excitations 
in a weakly interacting Bose gas (261 . The Hamiltonian 
for the combined system of an atom in a lattice site and 
the superfluid (for the motional atomic degrees of free- 
dom only) is given by 

H = //motion + //superfluid + //hit, (10) 

where //motion is a 3D harmonic oscillator Hamiltonian 
with frequency to, which describes the motional state of 
the atom, i/ S uperfluid is the Hamiltonian for the superfluid 
excitations, and H lnt is the interaction Hamiltonian. 

//superfluid = E + ^2 £ ( q ) kq^qi i 11 ) 
q#0 

where 6^ and 6 q are creation and annihilation operators 
for Bogoliubov excitations in the superfluid with momen- 
tum hq and energy e(q), and Eq is the ground state en- 
ergy of the superfluid. 

H^t = 9ab J 8p{r) <5p at om(r)d 3 r 

= g ab J Sp(r) S(r - f )d 3 r = g ab Sp(r), (12) 

where <5/3 a tom is the density operator for the motion of the 
atom, f is the position operator for the atomic motional 
states, Sp is the density fluctuation operator in the super- 
fluid, and g a b = 4TrK 2 a a b/(2p) is the coupling constant for 
the interaction, with a a b the scattering length for interac- 
tions between superfluid atoms and atoms in the lattice 
[22l | , and p = (m a mb) / (m a + rrib) the reduced mass of an 
atom in the lattice with mass m a and a superfluid atom 
with mass nib- The density fluctuation operator may be 
expressed as 8p = \&t\jr _ p where ^> — y/po + 8^ is the 
second quantised field operator for the superfluid and po 
is the mean condensate density. In terms of the creation 



4 



and annihilation operators for Bogoliubov excitations, we 
can write 



89 = -±= ( W q - r + v^e-^) , 



where V is the normalisation volume, 



1 



I -LI 



(13) 



(14) 



We define the projection operator, V , onto a basis di- 
agonal in the oscillator Fock states \m) {H mot - lon \m) — 
TiLo(m + 1/2)| m)) as 



VX = 



\m) (m| (m|X|m) 



so that 



Vw(t) = 2J \m)(m\pr, 



(19) 



(20) 



and 



L„ 



e q - (hq) 2 /(2m) — mu 2 



(15) 



The energy of excitations with momentum hq is 

s q = [u 2 (hq) 2 + (hq) 4 /(2m b ) 2 ] 1 ^, (16) 

and the speed of sound can be expressed as, u = 
V / .9h6Po/ m b; where gbb = ^irhabb / mb with abb the scatter- 
ing length for interactions between atoms in the super- 
fluid |26j . In a weakly interacting Bose gas at sufficiently 
low temperatures (where the condensate density is much 
smaller than the density of the normal component), the 
term from 59^ 89 may be neglected (see appendix lB)l . In 
this case we can write 



6p- 




+ « q )& q c iqr + K + v^e-^ 



For the motion of atoms in the lattice, we make the 
approximation that the damping of in each dimension 
can be can be considered independently, and thus we 
treat the atom as a ID oscillator with frequency u>, i.e., 



-ffatom = ^ I "'" + 2 



it; 



(18) 



where a is the lowering operator for the motional state 
of the at om. The p osition operator for the ID oscillator 
is x = \Jfij ' (2m a uj)(a + aJ), where m a is the mass of the 
atoms in the lattice. We can also write q.f — > q x x, where 
q x is the component of q in the direction of the oscillator 
motion. 



C. Damping Equations 

In deriving equations for the damping of the system we 
assume that the cooling rate is significantly slower than 
the period of the oscillator, and we treat the BEC as a 
reservoir in which the correlation time is much shorter 
than the correlation time of atoms in the lattice. Under 
these assumptions, the master equation for reduced den- 
sity operator describing the motional state of an atom 
in the lattice, w{t) = TrpJIU^)], where Tr^ denotes the 
trace over the superfluid states, is derived in appendixfXl 



Because we assumed that the oscillator trap frequency 
u> r _1 , where r is the characteristic timescale on which 
transitions take place due to interaction with the super- 
fluid, the coupling to off-diagonal elements of w(t) in the 
Fock state basis is very small, and the state occupation 
probabilities p n satisfy a closed set of equations. From 
appendix ^ we then see that 

Pm — ^ L n — >?nPn ^ F m — >n'Pm~ > r^^L[ n ^ rn (j) r i Prrl) • 
n>m n' <7n n 

(21) 

Here F n ^ m is the damping coefficient at zero tempera- 
ture for transitions from state n to state m, and H n>m are 
the coefficients of the finite temperature corrections due 
to the absorption and scattering of thermal excitations. 
The zero temperature damping coefficients are are 

F n ^ m = ~j-J2 l^n,m(q)l^(M" " m) - e q ), (22) 
q 

and the system matrix elements of the interaction Hamil- 
tonian are given by 



*(q) 



(m\H int \n) = ( Uq + Wq )*^( TO | e -^ £ | 



(23) 

Note that (1221) is Fermi's Golden rule for the transition 
from state |m) to state \n) via interaction with density 
fluctuations in the superfluid, and that the restrictions 
on the summation in l|21|) are written for clarity, but 
actually result from the delta function in (|22|l . due to 
which F n ^ m is only nonzero when n > m. 

The finite temperature corrections are given by 

H m , n = yY1 N{q)\Z n , m (q)\ 2 6(hLu\n - m\ - e q ), (24) 



where iV(p) = (exp[e p /(fcsT)] — 1) _1 is the mean num- 
ber of thermal Bogoliubov excitations with momentum 
hp present in the superfluid, with T the temperature in 
the normal component of the gas and fcs the Boltzmann 
constant. Because we are interested in cooling the sys- 
tem to its ground state, we assume that the temperature 
of the superfluid is small, ksT <C hoj. This is very re- 
alistic experimentally, and in this regime we can make 
the approximation H n m « 0. The consequences of these 
terms are described in section TlV CI 
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Also note that in the derivation of these equations we 
assume that the terms in Sp arising from 6^5^ may 
be neglected. This approximation holds when the tem- 
perature is much smaller than the critical temperature, 
so that the condensate density, po, is much larger than 
the density of the thermal component. The contribution 
from these terms at finite temperatures is estimated in 
appendix [5] 

From the transition rates, the energy dissipation rate 
of an oscillator in the state \n) can then be calculated as 

e(ri) = ^2?iu>(n-m)H m , n - ^ Koj(n-m)F n _> m . (25) 

m m<n 

The total energy dissipation rate for an atom in a mixed 
state can be written as e = Y] htonp n , which in terms 
of 123 is given by e = J2n z( n )Pn- 



D. Supersonic and Subsonic Motion Regimes 

We note that for typical experimental parameters in 
the lattice and the BEC, hu> 3> mbU 2 /2. For example, a 
Rubidium BEC with density po ~ 10 14 cm~ 3 and scat- 
tering length abb ~ lOOao, where ao is the Bohr Radius, 
has mbu 2 /(2h) — 2tt x 3.7 x 10 2 s _1 , whilst typically 
for an atom trapped in an optical lattice, to ~ 2tt x 10 5 
s . Thus, as the maximum velocity of the atom may 
be estimated as \J2Tluj / 'm a , we see that for a typical ex- 
perimental system and m a ~ mft, the atom velocities 
are supersonic. In this strongly supersonic regime the 
requirements of energy conservation in ill'l'li mean that 
even for a transition between states where m and n dif- 
fer only by 1, the excitations are in the particle branch 
of the Bogoliubov excitation spectrum. In this regime, 
the momentum of excitations generated in the superfluid 
hq ^s> mbU as h 2 q 2 /(2nib) > Uu) 3> rribii 2 /2. Hence, 



k q 2 /(2mb), and \u q +Vq\" 



1. 



If the superfluid was made sufficiently dense or strongly 
interacting, or the oscillator frequency u> was made suf- 
ficiently small that the motion of the oscillating atom 
was subsonic for all oscillator states which are initially 
excited, then energy conservation would cause the ex- 
citations to be in the phonon branch of the spectrum. 
In this regime, the momentum of excitations generated 
in the superfluid hq -C mt,u, so that e q ~ Tiuq, and 



+ v q \ 2 « hk/(2mbu). 



Note that the coefficients u q and v c 



can be related 



to the dynamic structure factor S(k, Co) of the superfluid 
which is often used in relevant literature |27|. In terms 
of the symbols used here, 5(k,cD) = \uk + Vk\ 2 - In the 
same way as previously discussed, S(k, Q) w 1 for large 
values of k 3> rribu/h, whilst for small k, 5(k, Co) oc k. 

In the following, we treat both the supersonic and sub- 
sonic regimes. As discussed previously, the supersonic 
regime is the more relevant of the two in current exper- 
iments. However, the subsonic regime could be specifi- 
cally engineered in experiments, and provides an inter- 



esting comparison in terms of the physics of the damping 
mechanism. 



IV. RESULTS 

The matrix elements in (|23H can be expressed in the 
position representation as 



ip*(x)tp m (x)dx, 



(26) 



where Iq — y/h/(m a co) is the oscillator length, q x is the 
component of q in the direction of the oscillator, ip n (x) = 
e -x /Wo)H n (x/l )/sJh2 n n\^ is the position wavefunc- 
tion for the state |n), and H n (x) is a Hermitc Polynomial. 
Using the identity dxe~^ x ~ v ^ H m (x)H n (x)dx — 
2 n i/7r ml y n ~ m L™~ m (— 2y 2 ), which assumes m < n, we 
can express the matrix elements (for m < n) as 



m! l 2 a 2 li 



l 2 n 2 

2 



l 2 n 2 
L a x 



so that 



9abP0 



poo ^2 

/ q 2 dqd(hLo(n — m) — Eg)- — \u q + v, 
Jo l oq 



2-nh n\ 

l q/V2 



loq/V? 



(27) 



This expression can be further analysed separately in the 
supersonic and subsonic motion regimes, where the re- 
sulting behaviour is remarkably different. 



A. Supersonic Case 

Applying to l|27|l the approximations given in section 
IIII Dl for the case of supersonic motion yields the expres- 



F„ 



sj (n—m)(m b /m a ) 

sj (n—m)(m b /m a ) 
9lbP0 m b ™ ! 

irh 3 l V2 n\ ' 



(28) 



The dimensionless function = 
Trh 3 l V2F n ^ m /(g 2 b p mb) is plotted in Fig. and 
shows the dependence of F n — > m on n and m. It is 
immediately clear that for all m < n the transition rate 
coefficient is significant. In fact, for all states \n), the 
transition rate directly to the ground state is of the same 
order as all other allowed transitions. This corresponds 
to the atomic motion generating a rich distribution 
of superfluid excitations, which is characteristic of the 
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regime where the motion of the atom is supersonic with 
respect to the velocity of sound in the superfluid (see 
section HVBl for a comparisonj. 

If we consider the energy dissipation rate (|25|) for a 
system in state \n) in the low temperature limit (fe^T <C 
hcu), then we see that the largest contribution comes from 
transitions directly to the ground state or first excited 
state (Fig. |3J). In addition, the complicated excitation 
spectrum results in a non-exponential energy damping 
law, i.e., the energy dissipation rate for a state n is not 
proportional to the energy of the state. We find instead 
for m a — nib — m (Fig. 0J) that 



i{n) 



9abP& 



,3/2 



[s(n)} 



3/2 



(29) 



where a = 0.3, and e(n) — Tiujn is the energy of state 
n measured with respect to the ground state. The 
total energy with respect to the ground state is then 
£ = Z) n £ ( n )Pn, so that i = -a Y, n l £ ( n )] 3/2 Pn ~ -ae 3/2 , 
provided that e 3 ! 2 ~ J2 n £ ( n ) 3 ^ 2 P n - The t' lme depen- 
dence of the total energy is then approximately given by 




FIG. 2: The value of F^ m = TTh 3 l V2F n ^ m /(g 2 ab p m b ), 
showing the coefficients of the transition rate from state n 
to state m in the case of supersonic motion, as computed 
numerically from equation l|28|l with m a = rrib- 



e(t) 



1 



-Va(i = 0)+ Set/2 



(30) 



where a = -g 2 ab pQm 3 / 2 uj 3 / 2 a/{-Kh 5/2 ^/2). 

The non-exponential damping law we obtain here can 
be understood in terms of a simple classical argument 
for a "foreign" atom moving uniformly through the su- 
perfluid at a supersonic velocity. If a a b is the scattering 
cross-section for the foreign atom interacting with the su- 
perfluid, then the average number of collisions per unit 
time is poCfabp/fn a , where p/m a is the velocity of the lat- 
tice atom propagating through the superfluid. The mo- 
mentum of the excitation generated in a collision is q oc p. 
Because the motion of the foreign atom is supersonic, 
the energy of the excitation is approximately q 2 /(2mt), 
and the energy dissipation rate i oc po&abP 3 /m 2 oc e 3 / 2 , 
which is the same energy dependence that we observe 
here. 

The slowest point in the cooling process is the cooling 
from the first excited state to the ground state, which 
is also the most important case for the low energy exci- 
tations which are likely to arise in quantum computing 
applications. Numerically we find that for m a = mb, 
F,' = 0.3789. Thus, 



v1.5< 



T 



I 

c 



'0.5 









-•- n=15 






-— n=10 






— n=5 




\ v 
V N, » 
— it % * 




><. X 








\ >v 








\ it 











10 



15 



m 



FIG. 3: The value of (n — m)i^_ >m , showing the contribu- 
tions to the energy dissipation of the system from transitions 
from state n to state m in the case of supersonic motion. 
These results are computed numerically from equation 
with m a ~ nib. 



Fi^o = 0.3789 



9lbPo m 
7rh 3 l V2' 



(31) 



The characteristic time for the transition from the first 
excited state to the ground state is then expressed in 
terms of the number of cycles by 



2tt 



1 



h 3 l uj 



0.3789 V2g 2 abPo m 

1 1 1 dab 

0.3789 16\/2 7r 2 p a 3 ab l ' 



(32) 



assuming that m a pa nib- In experiments, lo will typically 
be an order of magnitude larger than a a b, and the param- 



eter po«a6 



io- 



so the characteristic transition time 



from the first excited state to the ground state will be 
of the order of 10 cycles. It is interesting that the small 
prefactor in this expression is very important in giving 
such a rapid cooling rate. This rate is sufficiently fast 
to be useful experimentally, particularly given that the 
transition rates from states with higher quantum number 
to the ground state are all of the same order. 
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FIG. 4: The value of de'/dt = Y,t m<n ) ( n ~ m )K-> m plotted 
as a function of n for supersonic motion, showing the total 
rate of energy dissipation for a system instantaneously in a 
oscillator state with quantum number n. The points show the 
values computed numerically from equation 1281 with m a — 
nib, and the solid line is a fitted curve of the form de'/dt = 
an 3/2 , with a = 0.301. 



Subsonic Case 



Applying the approximations given in section Till Dl for 
the subsonic case to (|27|l . we obtain 



lQu{n— m) / (uy/2) 



lQU){n—m) /(uy/2) 



d£e-«Y (n - m) \L^ m (f)| 



(33) 



Fig.|5]shows F„_> m = ^■nm b hu b F n ^ rm /(g 2 ab p a Lj i ) plotted 
as a function of n and to. In contrast to the supersonic 
case, we see that i 7 ' n _ >m is very sensitive to the difference 
(n — to), and for sufficiently small hui/m a u 2 , the only 
significant contribution to transitions from the state |n) 
are transitions to state \n — 1). This can be seen very 
clearly in Fig. EJ which shows the contributions to the 
overall energy dissipation from the state |n). 

If we investigate the rate of energy loss e(n) from the 
state \n) as given by J2U, the only significant contribu- 
tion comes from the term where to = n — 1. For very 
small hui I m a u 2 , we can then expand the integrand near 
£ = 0, and, noting that L„_ 1 (Q) = n, obtain 



e{n) 



g 2 abPo u; 3 (n-l)\2l 2 L0 2 
,! 6u 2 



l2TTm a mbU 7 



hum 



-hojn. 



(34) 



Thus, as the energy of state |n) measured with respect to 
the ground state is s(n) = hum, the energy damping law 
is exponential. This is a direct consequence of the fact 
that, in contrast to the supersonic case, only the decay 



mode into the next lowest oscillator state is significant 
in the damping process. Fig. [7| illustrates from numerical 
calculations the linear dependence of the damping rate on 
n. Note that damping still occurs in this regime despite 
the fact that the velocity of the atom is slower than the 
speed of sound in the superfluid. This apparently contra- 
dicts the Landau derivation of the critical velocity in the 
superfluid. However, we note that the Landau criterion is 
a thermodynamic argument, and cannot be applied here, 
as the motion is accelerated. In fact, this damping law 
has an analogy with that for dipole radiation in classical 
electrodynamics (see section IVB|) . 

An analogy also exists between the supersonic and 
subsonic motion regimes here and regimes of large and 
small Lamb-Dicke parameter respectively in the context 
of laser cooling of trapped ions in a harmonic potential. 
In that system the Lamb-Dicke parameter, 77, is the ratio 
of the size of the ground state wavefunction to the wave- 
length of the cooling laser, and the interaction Hamilto- 
nian for the system is proportional to e lkx — e 171< ' a+a ^ 
Q. Here the interaction Hamiltonian is proportional to 
gigasa^ anc j wm i s t jy is a fixed parameter and q x is not, q x is 
constrained to be a small or large parameter by the con- 
servation of momentum when excitations are generated 
in the superfluid. 

During the cooling process, the coupling that ex- 
ists between two states, which is proportional to 
|(m| exp(— i77x)|n)| 2 is then analogous in the two cases. 
We observe cooling directly to the ground state from 
all excited states in the supersonic regime, and this 
is also a characteristic of cooling schemes in ion traps 
with a large Lamb-Dicke parameter. When r\ or q x are 
small (the subsonic regime or small Lamb-Dicke param- 
eter limit), the matrix elements simplify for m =/= n, 



(m I exp[i?7(a 



(m\ir)(a + a')\n), and coupling 



only exists between nearest neighbour states (this is 
known in ion trap cooling as coupling to the red and 
blue sidebands only). 



C. Finite Temperature Effects 

At finite temperatures, the terms proportional to H n ,m 
in (|21|) contribute heating effects due to the absorption of 
thermal excitations in the superfluid. When the temper- 
ature is significant, the final equilibrium motional state 
distribution will contain non-zero excited state probabil- 
ities. These can be calculated using the detailed balance 
condition [2^. Considering the transfer rates for atoms 
between oscillator states with consecutive quantum num- 
bers, we write for the equilibrium probability distribution 
Pn = p n (t -> 00), 



Fn+l^nPn+l — Hn+l,n(Pn ~ Pn+l) , 



so that 



H. 



Pn+l 



n+l,n 



(35) 



(36) 
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FIG. 5: The value of = 4Trm b hu 5 F n ^ m /(gl b p u} 3 ), 

showing the coefficients of the transition rate from state n to 
state m in the case of subsonic motion, as computed numeri- 
cally from equation 1331 for Iqui/u = 0.01. 
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FIG. 6: The value of (n showing the contribu- 

tions to the energy dissipation of the system from transi- 
tions from state n to state m in the case of subsonic motion. 
These results are computed numerically from equation [l.'i.'iti 
for Iqu/u — 0.01. 



Substituting the expressions from (1221) and l|24|l . and in- 
tegrating over the modulus of q, this expression simplifies 
to 



Pn+1 



N(qi) 
N( qi ) + 1 



Pn, 



(37) 



where fiqi — \/2hmhU) is momentum of excitations with 
energy e Ql = hcu. Thus, p n = [N ( Ql ) / {Nj qi ) + l)] n p , 
and using the normalisation condition Y^LoPn = 1, we 
obtain 



Pn = Pod 



-nhui/(k B T) 



1-e" 



huj/(k B T)\ -nhuj/(k B T) 




(38) 



FIG. 7: The value of de/dt — Y)/ m<n \ (n—m)F n ^, m plotted as 
a function of n for subsonic motion, showing the total rate of 
energy dissipation for a system instantaneously in a oscillator 
state with quantum number n. The points show the values 
computed numerically from equation 1331 . and the solid line 
is a fitted straight line of the form de/dt = an, with a = 
3.40 x 10" 5 . 



Hence, the equilibrium state occupation probabilities are 
simply given by the Boltzmann distribution, and the 
probability that an atom is in the ground motional state 
is p = l- e ~^/( feBT ). Provided k B T < Tico, the absorp- 
tion of thermal excitations will not significantly decrease 
the cooling rate, and will not prevent the cooling of es- 
sentially all of the population to the ground state. This 
obtainable under reasonable experimental conditions, for 
example, if lu ~ 2tt x 10 5 s -1 , hcu/kB ~ 5/jK, so that for 
T = 500nK, we then obtain 1 - p w 5 x 1CT 5 . 



V. THE SEMI-CLASSICAL APPROXIMATION 

A. Supersonic Case 

It is interesting to compare the fully quantum cal- 
culation of the damping rates to the calculation in the 
semi-classical approximation. Using this approximation, 
the calculation is performed similarly to the calculation 
of damping due to radiation from an oscillating charge, 
which provides a useful physical analogy between the two 
situations. 

In this calculation we make use of the relationship be- 
tween quantum matrix elements and the Fourier com- 
ponents of the classical trajectory of the system [2^ |. 
Strictly speaking, this approximation is valid only when 
the equivalent quantum matrix elements are taken be- 
tween states of large quantum number, and where the 
difference in the quantum numbers is small relative to 
the quantum numbers. We will discuss the validity of the 
approximation in practice at the end of the calculation. 
The classical trajectory of the atom in the lattice may 
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be written in ID as r(t) — > r max cos(wi)z, where z is the 
axial unit vector along the lattice. Because the motion 
is periodic with period 2w/ui, the frequency spectrum of 
the resulting excitations will be discrete with frequencies 
urn for integer n. Analogously to (|22|l . we then write the 
rate of energy dissipation for the atom in the lattice (at 
zero temperature) as 



e = — — \Tq(um)\ 2 5(hujn — e q )huin, 



(39) 



where 



^ |T q (wn)| 2 = X 



q 



N f 



2k/u 



(40) 

with \Nf) the final state of the superfluid (normally a 
state with a particular number of excitations of momen- 
tum hq) . This expression is also averaged over the initial 
state of the system |iVj), which will usually correspond 
to a thermal distribution of excitations. 

Assuming that we are in the supersonic motion regime 
and applying the approximations given in section fill Dl 
we obtain 



T q (ujn) 



9abJP0 w 



2tv/u 



W 27T Jlt 

Using the identity 



-ifefmaa cos(t)g— iumt ^j. 

(41) 



1 

27 



2n 



-\z cos(C) p -inC 



dC 



Jl{z), (42) 



where J n (z) is an ordinary Bessel Function, and integrat- 
ing over the angular values of q in spherical coordinates 
then gives 



9 2 ab Po 



]T / dqq 2 / d£J 2 n 



2tt 

xd(hujn — e q )Lon. 
We now integrate over q to give 



(43) 



2 3/2 3/2 

9abP0 m b <± ' 



£n 3 / 2 f d^(WH), (44) 



where a = r max y/2mbU) /%. We can see that many values 
of n contribute significantly to this sum, which is analo- 
gous to the full quantum result, in which many different 
transitions between oscillator levels had significant coef- 
ficients F n ^ m . As noted in section llV Al this fact arises 
from the the motion of the oscillating atom being faster 
than the speed of sound in the superfluid. This spectrum 
of generated excitations can be seen as being analogous 
to the result for electromagnetic radiation from a charge 
moving faster than the speed of light (in a dielectric), 
which can be computed semi-classically using a similar 
method to that used here. 



It is possible to determine analytically the functional 
dependence of l|44|) on r max by finding an approximate 
expression for the integral over £. In the limit where the 
argument of the Bcsscl function is large, we can write 



F(a,ri) 



d£J 2 (fry/*) 



2 f 1 dj 2 cos2 (W" - - tt/4) 
J to 



Co 6 Tovn 



where £o = \/n/a is the lower limit for £ in which 
the cosine approximation of the Bessel function is valid. 
This expression is strictly only valid for n -C a 2 = 
2r 2 nax mf ) u;/?i. At larger values of n, F{a,n) is ex- 
ponentially small, and the functional dependence of 
n 3 / 2 F{a, n) on a can be found from the point at 
which the summation is cut off, and for a system of en- 



ergy WaWfL/ 2 . «n 



m a wr 1 2 nax /(2?i) = a 2 m /(4mi). 



Approximating the sum by an integral, we can then write 



-C 



g 2 b p m 2 mb[l + 21n(4m 6 /TO a )] 



,3^,3 



327r 2 n 4 



(46) 



where C is a constant which for large values of a is in- 
dependent of a. Fig. [S] shows a numerical calculation of 
C(a), from which we observe that for large a, C ~ 1.75. 
Moreover, the approximation is also very good for small 
values of a > 2, so that C is essentially a constant for all 
physical values of a. 




FIG. 8: The value of C computed numerically as a function 
of a = \/2?"max/'o by comparison of the results from l|44|l 
and 146H . Note that this curve is discontinuous because of 
the discrete sum in jUt i which was cut off at the highest 
integer less than a, and that C = for a < 2, because a < 2 
corresponds to a sum cut off at n = 0. The value of this 
function in the limit as a — > oo gives C ~ 1.75. 



If we use the classical expression r max = ^/2s/(m a U} 2 ), 
where e is the energy of the oscillating atom, we can 
rewrite (I46|l as 

s = -C[l + 2H^b/ ma )] ^ 9 ^S mb ^' 2 - ^ 
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As in the quantum case, the damping is non-exponential 
as a result of the rich distribution of generated exci- 
tations and instead e oc e 3 / 2 . If we compare this re- 
sult to that from equation l|29[) . the ratio of the semi- 
classical result to the quantum result for m a = mi is 
C[l + 41n(2)]/(8a7r) w 0.88. The reason for this be- 
comes clear when we examine the terms of the series 
~^2 k k 3 / 2 F(2n, k) (Noting that if we begin in the initial 
state |n) then a = 2n), and compare them to the equiv- 
alent terms in the quantum calculation, kF n _ 
This is shown for an initial state n = 10 in Fig.El We see 
that the terms agree well for small k but that they di- 
verge as k — > n. This is because the equivalence between 
the semi-classical result from the Fourier spectrum and 
the quantum matrix elements is strictly only valid when 
k is small. Because in the calculation of energy dissipa- 
tion rates the terms are weighted by an additional factor 
of k, the terms where the largest discrepancy arises are 
always significant in the calculation of the damping rates, 
and thus this discrepancy does not significantly decrease 
asn-> oo. 
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FIG. 9: Numerical values of the quantum result -f\o^(io-fc) 
(solid line) and the semi-classical result F(20, k) (dotted line) 
in the supersonic regime. Note that we observe very good 
agreement for small k, but the results diverge for higher values 
of k. 



B. Subsonic Motion 

In addition to the approximations given in section 
1111 Dl we note that for the purposes of the semi-classical 
calculation in the subsonic regime, |q.r max | < qv meiK /u> = 

Vmax/U < 1- Thus, 



1 

2^ 



271 



-iq*r m „ cos Cg-inC ^ 



glbpo 



xS(fllU — £ q )TlU) 

Ylitu' 'mf,m a 



dqq 3 / d£ 



(49) 



2^1 q * r » 



: cos(C)e-^dC = i^^5 n , ± {,48) 



As mentioned in section IIVBI damping occurs here de- 
spite the fact that the velocity of the atom being slower 
than the speed of sound in the superfluid appears to con- 
tradict the Landau derivation of the critical velocity in 
the superfluid, and we obtain an exponential damping 
law. In the same sense that the previously discussed case 
of supersonic motion is analogous to radiation from a 
charge moving faster than the speed of light in a dielec- 
tric, this case is analogous to dipole radiation from an 
accelerating charge. The approximation made that re- 
sults in only one term in the sum being significant, 14811 . 
similarly corresponds to the dipole approximation in non- 
relativistic quantum electrodynamics. 

Note that if we substitute e = hcun into (QSJ, then 
we obtain exactly the same result we obtained from the 
quantum case (|34fl . The semi-classical approximation 
works extremely well here, because the only significant 
contribution to the quantum calculation comes from ma- 
trix elements between states with quantum numbers dif- 
fering by one. 



VI. IMMERSION IN A STRONGLY 
CORRELATED ID SUPERFLUID 

In this section we investigate the damping that occurs 
when the lattice is immersed in a quasi-one dimensional 
superfluid, which is an example of strongly correlated 
quantum liquid. In a real experiment this setup is not 
particularly practical for cooling the motion of the atoms. 
For a gas to be quasi-one dimensional, the excitation 
modes in the transverse directions must have energies 
larger than all other significant energy scales in the sys- 
tem, and so the oscillator energies for lattice atoms, tiu> 
must be much smaller than the energies of the transverse 
excitations in the superfluid. Furthermore, the motion 
of the oscillator will only be damped in one dimension 
(along the direction of the quasi-lD superfluid), and so 
the oscillator should be made strongly anisotropic so that 
in the transverse directions the oscillator is always in the 
motional ground state and need not be cooled. However, 
the study of the cooling process in this context is still 
interesting because the lattice atom in this setup could 
be used as a probe to provide spectroscopic information 
about the ID Bose gas. 

In general the excitation spectrum of such a one- 
dimensional Bose gas is complicated. In the case of short- 
range interactions between the particles exact analyti- 
cal solution exists both for the ground state wavefunc- 
tion and for the excitation spectrum |30j | for arbitrary 
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strength of the interparticle interactions and the excita- 
tion energies. However, in the limit of long wavelength 
the excitations are phonons and the system can be de- 
scribed within a hydrodynamic approach. Following |3l| 
we represent the field (Bose-particle annihilation) oper- 
ator in the form: ^>(x) cx vVo + Spe 1 ^, where <f> and 8p 
are phase and density fluctuation fields respectively and 
obey the commutation relation [Sp(x), 4>{y)} = \5(x — y), 
and Pq is the ID density (averaged, in practice, over the 
transverse directions). The low-energy effective Hamil- 
tonian for the liquid is then 



% 

2^ 



dx[vj{d x 0) 2 +v n {tt6p) 2 }, 



(50) 



where vj = irhpo/rrib, vn — k/ (irhpo), and k is compress- 
ibility per unit length. The excitation spectrum corre- 
sponding to this Hamiltonian satisfies a linear dispersion 
relation e q = hv s q, where the velocity of sound is given 
by v s = (vjVn) 1 / 2 . 

The parameters vj and vn are phenomenological and 
can be found from the exact Lieb-Liniger solution |3Cj . 
The dependence on the interaction strength between gas 
particles can be described using the dimensionless param- 
eter, 7 = mb9bb / (ft 2 Po) ■ In the week interaction limit, 
7 <C 1, the velocity of sou nd is given by the usual Gross- 
Pitaevskii value: v s — \J gbbPo/ m b- If the interaction 
is very strong, 7 » 1, then the interaction effectively 
makes the particles impenetrable, and hence in a true ID 
system, indistinguishable from Fermions. This is called 
the Tonks gas regime, and the sound velocity is equal to 
the effective Fermi velocity: v s — irhpo/mb. The energy 
spectrum is linear for e q <C gbbPo (the chemical potential 
of a weakly interacting Bose-gas) and e q <C nh 2 ' p\i '(2mf,) 
(the Fermi energy of the Tonks gas) for the cases of weak 
and strong interactions, respectively. At higher energies 
the excitation spectrum is no longer universal and de- 
pends on the details of the interparticle interactions. Be- 
cause e q and the trapping frequency u> in the lattice are 
related via energy conservation, the motion of the lattice 
atoms must then be subsonic with respect to v s for the 
model to be valid. 

The operator for density fluctuations in this regime is 
given by 



1/2 



Sje" 



iqx 



(51) 



where L is the length of the BEC and K = {vj/v^) 1 / 2 . 
The quantity K depends on the interparticle interactions 
and is related to the scaling dimension of the particle 
field operator: {&(x)i>{x')) ~ \x - x'^ 1 ^ 2 ^ for large 
\x — x'\. The function if (7) monotonically decreases as 



tt[7- (1/2tt) 7 3 / 2 ]- 1 /2 and 
Note also that for the quasi 



7 grows, so that if (7 — » 0) 
if(7-^oo)«(l + 2/ 7 ) 2 [3L 

one-dimensional system, gbb = 'inh 2 abb / 'rribl 2 ^, where lj_ 
is the transverse confinement length of the BEC, provided 
that a B < lj_ HI. 



In the limit of small oscillation frequencies lo, we apply 
the same approximation l|48|l used in section IV Bl and 
obtain 



dqq 



qr„ 



S(huj — hv s q) 



-alb 



Klu 2 



(52) 



For small 7, K 
^gbbPo/mb, so 



7r^V / / 9 o/( m b.9bb) and v s 



"9ab^ 



7/4 



fc 7/4 9/4 



e. 



(53) 



The transition rate constant is then r e ~ 
u(gab/gbb) 2 (hu/pogbb){mbgbb/h 2 po) 3/i (m b /m a )/y/ir < 
lo and hence is generally small. In the opposite limiting 
case for large 7, K w 1 and v s = irhpo/mb, so 



9 9 4 
~9ab U ' m b 

7r 5 h 5 m a pQ 



e. 



(54) 



Here, T e - uj(m b g a b/h 2 ' p ) 2 (ujmb/hpl)(m b /rri a ) /tt 5 . 
Thus, in this regime, the damping rates can be made 
very fast, provided that 7 = mbg a b/Ti 2 po is made very 
large. However, this regime is difficult to obtain experi- 
mentally, and in most current experiments 7 ~ 1. 

In both cases the damping that we obtain is exponen- 
tial, which again arises because the motion we consider 
is subsonic, and produces excitations at only one signif- 
icant momentum. The energy exchange rate grows as a 
function of uj, in a manner analogous to dipole radiation 
in quantum electrodynamics. 



VII. SUMMARY 

We have shown that the immersion of a system of 
atoms in an optical lattice in a superfluid causes damping 
of atoms in excited motional states, and that this damp- 
ing can be used to transfer these atoms to the ground mo- 
tional state whilst preserving their initial internal state 
and any entanglement between the atoms. For typical 
experimental parameters, this transfer occurs in a char- 
acteristic time of around 10 oscillator cycles, which is 
sufficiently rapid to be useful experimentally. These typ- 
ical parameters come from a regime in which the atoms 
in the lattice are moving faster than the velocity of sound 
in the superfluid, which generates a rich distribution of 
excitations, involving significant transitions from all lev- 
els directly to the ground state. In the opposite regime, 
where the velocity of the atoms in the lattice is signifi- 
cantly slower than the speed of sound in the superfluid, 
damping still occurs because the motion is accelerated, 
but only transitions between neighbouring oscillator lev- 
els contribute significantly to the damping process. 
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Provided that the temperature in the non-superfluid 
fraction of the gas is much smaller than the oscillator 
level spacing in the lattice, heating effects due to ab- 
sorption of thermal excitations is not a significant effect 
in this process. This is the case for experimentally re- 
alisable conditions. At higher temperatures, the system 
would be cooled not to the ground state, but to a thermal 
distribution of motional states corresponding to a Boltz- 
mann distribution with the same temperature as that in 
the normal component. 

The supersonic motion regime discussed here is readily 
realisable in present experiments. Together with a careful 
choice of internal atomic states used to encode a qubit, 
this damping mechanism thus provides a decoherence- 
free means to cool an atomic qubit to its motional ground 
state. 



Acknowledgments 

This work was supported in part by the Austrian Sci- 
ence Foundation FWF, E.U. Networks and the Institute 
for Quantum Information. 



APPENDIX A: DERIVATION OF THE MASTER 
EQUATION 

We treat the superfluid with Bogoliubov excitations as 
a reservoir, with density operator R. In the interaction 
picture, and after making the Born-Markov approxima- 
tion, the master equation for the density operator w of a 
system which interacts with a reservoir via an interaction 
Hamiltonian Hi n t can be shown to be given by |33| 

1 r* . 

w = --^ dt'Tr R [H int (t),[H iat {t'),w(t)®R]}, (Al) 
n Jo 

I 



where Tr R denotes the trace over the reservoir states. 

We write s q ,i = e 19x2; , s qi2 = e -19 ^, f q ,i = 6 q and 
r qj 2 = bt, so that 

q i=l,2 

and then substitute this expression into QAlj l to give 




[Si{t)Sj(l/)w{t) - Sjit^W^Siit)} (f<(t)f jtf)) H 

+ [«;(*)«,•(*%(*) - Si(t)«j(t)« i (*')] {ttfftiffiU, 

(A3) 

where we have used the cyclic property of the 
trace, dropped the operator subscript q, and written 
Ttr(RA) = (A) R . We have also used the fact that 
(f q ,i(Of q ;,i(t))ii = for q^q'. 

Proceeding in the standard way, we change the vari- 
able of integration to r = t — t' , and note that 6 q (i — 
t) = C - i6bT / n b cl (t)e i " bT / n = e i£ i T / ri 6 q (t) and similarly 
g-i^r/figt (£) e ii W fi = e^^H^t). 

We then make use of the assumption that w(t) as 
Vw{t) (see H19(l 'l. and write the master equation in a Fock 
state representation. Noting also that (bq(t)bq(t)) R — 
(S q (f)S q (t)) r = 0, we obtain 



w 



2 9lbP0( U q + V qY 

vn 2 



X/X/ / dr l m )( m l e iq!c:i \n)(n\e iq * i: \m)(m\p m e 



\Ljr{m—n) 



+ \m)( m \e- i ^\n}(n\e- i ^- £ \m)(m\p n e i ^ n ~ m ^ ' h {bjj\) R + c^ T ' h {b%) R ) . (A4) 

Assuming that the correlation time of the superfluid reservoir is much shorter than that in the system we can 
extend the integration over r — > oo, and making the replacement j °° dre l( - £_eo -' r / 7i — > nhS(e — eq), we obtain 



27r gl b Po(u q + v q ) 

vn 



y^y^ \{m\e iq * x \n)\ 2 j[<5 {%u{n - m) - e q )p„ - 5 {%w{m - n) - e q )p m ] (& q 6 q )ii 



+ [5 (hu(m - n) - eq)p n - S {huo(n - m) - e q )p m ] (6 q 6 q )fl| 



(A5) 



The first two terms here (those proportional to (& q & q )fl) superfluid, whilst the second two terms (those propor- 
describe the damping by creation of excitations in the tional to (6 q 6 q ) R ) describe heating effects by absorption 
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of thermally generated excitations. At finite tempera- 
tures, the reservoir correlation functions are given by the 
number of thermal excitations iV(q) with momentum Tiq, 
(b&J a = N(q). 



APPENDIX B: ESTIMATION OF S&Sty TERMS 

The heating effects due to absorption of thermal exci- 
tations has already been discussed in section llV CI where 
the equilibrium distribution at finite temperatures was 
shown to be a Boltzmann distribution. The small addi- 
tional damping terms arising at finite temperatures from 
the S^S^> term, which are small when the condensate 
density, po is large and which were omitted when the den- 
sity fluctuation operator Sp was originally written, may 
be estimated using a semi-classical treatment. The oper- 
ator for the additional density fluctuation terms is given 
by 

1 



Sp' = = TfYl M P u P'«P«p' ei(p ~ p,) ' r 

p,p' 



+ Vr>V 



/fit a D /e- i < p - p '>- r + ■ 



p Up' IX p U p 

+ u P vap4' e_i(P+P ' ) ' 



u p Vp> a p a p ie 



(p+p') 



(Bl) 



The first two terms in this expression correspond to the 
inelastic scattering of thermal excitations with momen- 
tum hp to excitations with momentum hp' , and the sec- 
ond two correspond to the absorption and emission re- 
spectively of two excitations with momenta hp and hp' . 

For the case of supersonic motion where u q — > 1 and 
v q — ► 0, the correction to the dissipation rate is then 
given by 



££[tf(p)-JV(p')]*(& 



ton - 



p,p' n 

/>27T I ' UJ 

^ / e i(p-p').r(t) df 

2W 



hum, 



£p> + S p ) 
(B2) 



where, as before, iV(p) = (exp[e p / (ksT)] — 1) _1 is the 
mean number of thermal Bogoliubov excitations with 
momentum hp present in the superfluid. 

In order to cool the system to the ground state we al- 
ready require %u> 3> fcgT, which has been shown to be 
a reasonable experimental condition in section TlV CI In 
this case, the thermally generated excitations with mo- 
mentum hp will have a much smaller energy than the 
scattered excitations, which have momentum hp' . Also, 
e p / > huj > k B T and AT(p') w 0. Thus, 



LOU — E p i ) 



P,p' n 

2ir 



hun, 



(B3) 



This result is proportional to the density of thermal ex- 
citations and essentially describes the classical friction 
due to scattering of thermal excitations by the moving 
particle. 

If ksT -C to;,m 2 /2, then e p « hup. The additional 
damping is then given in terms of the rate e in (|47|) by 



eC(3) 



2n 2 po(hu) : 



-{k B Tf 



(B4) 



where Q{x) denotes the Riemann Zeta function. Note 
that because the wavenumber of phonons in this regime 
is of the order of k B T / '(hu), this result is proportional 
to the density of thermal phonons, Pphonons- Thus the 
additional damping term is equal to that in l|47l) . but 
with the numerical coefficient modified, and the density 
of the condensate po replaced by the density of thermal 
phonons, p p honons- This term will always be small, as in 
this regime T < T c , the critical temperature of the Bose 

gas, SO p Q > /Ophonons- 

If k B T > m b u 2 /2, then e p 
additional damping is then 



h 2 p 2 /(2m b ). The rate of 



, _ i((3/2)m 3 b /2 3/2 

£ wh^p-tf {kBT) ■ 



(B5) 



For a uniform Bose gas the critical temperature for Bose 
condensation can be expressed as |3 II 



kuT r 



2nh 2 p 2 t /3 
m b [C(3/2)]2/3' 



(B6) 



where pt = Pa + Pn is the total density, and p n is the 
density of the normal component, so that we can rewrite 
(lB5l) as 



ept 
2 Po 



3/2 



(B7) 



where we have used the well known result p n = 
p t (T/T c ) 3 / 2 j^ll- Thus, this result has the same form as 
the damping rate obtained in (|47|) . but the condensate 
density is replaced by the density of the normal com- 
ponent, and the numerical coefficient is decreased by a 
factor of 2. Again, at small temperatures compared with 
the critical temperature, T <^ T c , when p n <C pa, the 
contribution from this term will be small. 

The same calculation can be performed for the sub- 
sonic case. In this regime, the contribution from the 
terms involving a p a p r and Spap» is small, because the 
double summation over p and p' is restricted by energy 
conservation such that |e p + gy| = hton, and in the sub- 
sonic case, this quantity is always small. With respect to 
the subsonic energy dissipation rate in l|49() . i, we obtain 



480 i oom b u 5 ?i' 



-\k B T)\ 



(B8) 
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Note that as Tiuin <C mi,u 2 /2, this expression is derived 
considering only the case where kgT <C m^v? /2. It can 
be shown that in the limit ksT <C mu 2 /2 that the density 
of the normal component p n is given by |3J| 



Pn 



so that we can write 



2n 2 {k B T) A 
45m b h 3 u 5 

8J as 



(B9) 



Again, this result is a modification of the zero- 
temperature damping result, with the condensate density 
replaced by the density of the normal component and the 
numerical coefficient decreased. In the limit T <C T c , as 
with the supersonic results, this result will be small, as 
Pn < Po- 



3£p n 
64p ' 
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